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I Abstract 

(N ; 

^ ■ We systematically analyze constraints on supersymmetric theories imposed by 

^ . the experimental bounds on the electron, neutron, and mercury electric dipole 

^ I moments. We critically reappraise the known mechanisms to suppress the 

I EDMs and conclude that only the scenarios with approximate CP-symmetry 

^ ' or fiavour-off-diaeonal CP violation remain attractive after the addition of the 

O ■ 

■ mercury EDM constraint. 

^' 

^ ; 1 Introduction 

> ■ 

• T-j , There are a number of reasons to suspect that there are additional sources of CP viola- 

I tion beyond those of the Standard Model (given by 9 and 6km)- The most compelling 

a I one is that the SM is unable to explain the cosmological baryon asymmetry of our uni- 

verse. Also, the Standard Model is very unlikely to be the "ultimate" theory of nature. 
Most extensions of the SM bring in new sources of CP-violation. In particular, the 
most attractive one - the Minimal Supersymmetric Standard Model (MSSM) allows 
for new sources of CP violation in both supersymmetry-breaking and supersymmetry- 
conserving sectors, and there are no compelling arguments for them to be zero. In 
addition, the improving precision in the measurements of the CP-observables such as 
Acp{B — i> ipKs) [in may soon reveal deviations from the SM predictions. 

The most stringent constraints on models with additional sources of CP violation 
come from continued efforts to measure the electric dipole moments (EDM) of the 
neutron 0, electron 0, and mercury atom |^ 

dn < 6.3 X 10"^^ e cm (90%CL), 
4 < 4.3 X 10"^^ e cm , 
dHg < 2.1 X 10"^^ e cm . (1) 



1 



With the expected improvements in experimental precision, the EDM is hkely to be 
one of the most important tests for physics beyond the Standard Model for some time 
to come, and EDMs will remain a difficult hurdle for supersymmetric theories if they 
are to allow sufficient baryogenesis. Indeed it is remarkable that the SM contribution to 
the EDM of the neutron is of order 10~^° e cm, whereas the "generic" supersymmetric 
value is 10~^^e cm. 

In this paper we analyze neutron, electron, and mercury EDMs in the context of 
R-parity conserving supersymmetric theories. In particular, we reconsider the known 
mechanisms to suppress EDMs in light of the recently reported bound on the mercury 
EDM 1^. These include SUSY models with small CP phases, models with heavy 
sfermions, the cancellation scenario, and models with flavour off-diagonal CP violation. 
We also study to what extent different scenarios rely on assumptions about the neutron 
structure, i.e. chiral quark model vs parton model. 

The paper is organized as follows. In section 2 we present general formulae for the 
EDMs and discuss their model-dependence. In section 3 we define our supersymmetric 
framework and present all relevant supersymmetric contributions to the EDMs. In 
particular, we analyze the importance of the two-loop Barr-Zee and Weinberg type 
EDM contributions. Section 4 is devoted to the study of the EDM suppression mecha- 
nisms. First we consider in detail the "canonical" scenarios: suppression due to small 
CP phases or heavy sfermions. Second, in the context of the cancellation scenario, we 
analyze the possibility of the EDM cancellations in two classes of models: mSUGRA- 
like models with nontrivial gaugino phases and D-brane models. Finally, we discuss 
the EDM suppression in models with flavour-off- diagonal CP violation. In addition, 
we present new model-independent bounds on the sfermion mass insertions imposed by 
the electron, neutron, and mercury electric dipole moments. In section 5 we overview 
and discuss our results. 

2 Electron, neutron, and mercury electric dipole 
moments 

Let us flrst summarize the contributions to the three most signiflcant EDMs, beginning 
with the most reliable, the electron EDM. 

2.1 Electron EDM 

The electron EDM is deflned by the effective CP-violating interaction 



where F^'^ is the electromagnetic fleld strength. The experimental bound on the elec- 
tron EDM is derived from the electric dipole moment of the thallium atom and is given 



(2) 



by i 



de < 4 X 10 e cm . 



(3) 
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In supersymmetric models, the electron EDM arises due to CP-violating 1-loop dia- 
grams with the chargino and neutralino exchange: 

4 = df + df . (4) 

Since the EEDM calculation involves little uncertainty it allows to extract reliable 
bounds on the CP-violating SUSY phases. 



2.2 Neutron EDM 

The neutron EDM has contributions from a number of CP-violating operators involv- 
ing quarks, gluons, and photons. The most important ones include the electric and 
chromoelectric dipole operators, and the Weinberg three-gluon operator: 

— -^d^ f abcG a^ipG^j^G c\cr^^^^" i (5) 

where Ga^v is the gluon field strength, T°- and fahc are the SU(3) generators and group 
structure coefficients, respectively. Given these operators, it is however a nontrivial task 
to evaluate the neutron EDM since assumptions about the neutron internal structure 
are necessary. In what follows we will study two models, namely the quark chiral model 
and the quark parton model. Neither of these models is sufficiently reliable by itself , 
however a power of the combined analysis should provide an insight into implications of 
the bound on the neutron EDM and in particular comparing them gives some indication 
of the importance of these systematic errors in, for example, cancellations. A better 
justified approach to the neutron EDM based on the QCD sum rules has appeared in 
P] and earlier work P|, |jl3|. We note that in any case the NEDM calculations involve 



uncertain hadronic parameters such as the quark masses and thus these calculations 
have a status of estimates. The major conclusions of the present work are independent 
of the specifics of the neutron model. 

i. Chiral quark model. This is a nonrelativistic model which relates the neutron 
EDM to the EDMs of the valence quarks with the help of the SU(6) coefficients: 

4 1 

dn = -^dd - -du . (6) 
The quark EDMs can be estimated via Naive Dimensional Analysis [H as 



d, = r]'^df + + V^'^d'' , (7) 



where the QCD correction factors are given by rj^ = 1.53, rj'^' c:^ rj'-^ c:^ 3.4, and 
A ~ 1.19 GeV is the chiral symmetry breaking scale. We use the numerical values for 
these coefficients as given in R]. The parameters t]'^''^ involve considerable uncertainties 
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steming from the fact that the strong couphng constant at low energies is unknown. 
Another weak side of the model is that it neglects the sea quark contributions which 
play an important role in the nucleon spin structure. 

The supersymmetric contributions to the dipole moments of the individual quarks 
result from the 1-loop gluino, chargino, neutralino exchange diagrams 

^E,C ^ (E,C) ^ (E,C) ^ (E,C) ^ (g) 

and from the 2-loop gluino-quark-squark diagrams which generate d'-'. 



ii. Parton quark model. This model is based on the isospin symmetry and 
known contributions of different quarks to the spin of the proton [|1^]. The quantities 
Aq defined as (?T-||^7/^755'|n) = 5*^, where 5*^ is the neutron spin, are related by 
the isospin symmetry to the quantities (Ag)p which are measured in the deep inelastic 
scattering (and other) experiments, i.e. = {Ad)p, A^ = (Au)p, and A^ = {As)p. To 
be exact, the neutron EDM depends on the (yet unknown) tensor charges rather than 
these axial charges. The main assumption of the model is that the quark contributions 
to the NEDM are weighted by the same factors Aj, i.e. |]10| 



d„ = r/^(A,df + AX + A.rff). (9) 
In our numerical analysis we use the following values for these quantities A^^ = 0.746, 



A„ = —0.508, and As = —0.226 as they appear in the analysis of Ref . ||1 1|| . As before, 
we have 

< = 4(^) + rff + (10) 

The major difference from the chiral quark model is a large strange quark contribution 
(which is likely to be an overestimate 0). In particular, due to the large strange 
and charm quark masses, the strange quark contribution dominates in most regions 
of the parameter space. This leads to considerable numerical differences between the 
predictions of the two models. 

The current experimental limit on the neutron EDM is 

dn < 6.3 X 10~2^e cm . (11) 

2.3 Mercury EDM 

The EDM of the mercury atom results mostly from T-odd nuclear forces in the mercury 



nucleus |T2[, which induce the effective interaction of the type (I ■ V)5(r) between the 
electron and the nucleus of spin I In turn, the T-odd nuclear forces arise due to the 
effective 4-fermion interaction ppni'^^^n. It has been argued ^ that the mercury EDM 
is primarily sensitive to the chromoelectric dipole moments of the quarks and the limit 

i 

dng < 2.1 X lO^^^e cm (12) 

can be translated into 

- rf^ - 0.012rff 1/^, < 7 X 10-2^ , (13) 
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where Qs is the strong couphng constant. As in the parton neutron model, there 
is a considerable strange quark contribution. The relative coefficients of the quark 
contributions in (0) are known better than those for the neutron, however the overall 
normalization is still not free of uncertainties [|T3[| . 



3 The EDMs in SUSY models 

We will study supersymmetric models with the following high energy scale soft breaking 
potential 

- A^Yi; H2qLiU*Rj + h.c.) + ^{msgg + m^Ww'^ + m,BB) , (14) 

where denotes all the scalars of the theory. We generally allow for Ai ^ ^ Ad 
as well as nonuniversal gaugino and scalar masses, which is important for the analysis 
of D-brane models. The /x, 5, A^, and mj parameters can be complex, however two 
of their phases can be eliminated by the U{1)r and U{1)pq transformations under 
which these parameters behave as spurious. The Peccei-Quinn transformation acts 
on the Higgs doublets and the "right-handed" superfields in such a way that all the 
interactions but those which mix the two doublets are invariant. The Peccei-Quinn 
charges are Qpgifi) = QpgiBfi), Qpq{A) = Qpgirrii) = 0. The U{1)r transforms 
the Grassmann variable 6 — >■ 6e^°' and the fields in such a way that the integral of the 
superpotential over the Grassmann variables is invariant, i.e. the f/(l)_R charge of the 
superpotential is 2. As a result, QR^Bfi) = Qr{^i) — 2, Qr{A) = QR{mi) = —2. The 
six physical CP-phases of the theory are invariant under both U{1)r and U{1)pq, and 
can be chosen as 

Arg{A*imi) , Arg{{Bfi)*fiA^) , (15) 

where i = 1, 2, 3 and a = d,u, I. All other CP-phases can be expressed as their 
linear combinations. If the A-terms are universal, there are four physical phases 
Arg{A*mi) , Arg{{Bfi)*fiA). 

It is customary to choose the phase convention in which the Higgs potential param- 
eter Bn is real. In this case, the physical phases become Arg{A'^mi) and Arg{fiAa). If 
universality is assumed, the number of physical phases reduces to two. In what follows 
we will set m2 to be real by a U{1)r rotation. 

Nonuniversality will play a crucial role in the D-brane and flavour models' analysis, 
but otherwise does not lead to different conclusions for the models we study. Thus we 
will assume universal A-terms and gaugino masses unless otherwise specified. 

In what follows we use tan/3, mo. A, rrti as input parameters and obtain low energy 
quantities via the MSSM renormalization group equations (RGE). We also assume 
radiative electroweak symmetry breaking, i.e. that the magnitude of the /i parameter 
is given (at tree level) by 



1^1 = tan^/3-1 2"^^ • 
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The phase of /i is an input parameter and is RG-invariant. Our numerical results are 
sensitive to the quark masses which we fix at the Mz scale to be: m„. = (0.005, 1.40, 165) 
GeV and m^i. = (0.010, 0.194, 3.54) GeV. The light quark masses are poorly determined 
(in fact iTLu = is not excluded) which results in the uncertainties of the EDM nor- 
malization; for definiteness, we have chosen the light quark masses as they appear in 
0. The GUT scale is assumed to be 2 x lO^^ GeV. 

It is well known that 0{1) supersymmetric CP phases generally lead to unaccept- 
ably large electric dipole moments which constitutes the SUSY CP problem. In this 
paper we consider different mechanisms for suppressing EDMs and analyze them in 
detail. 



3.1 Leading SUSY contributions to the EDMs 

In this subsection we list formulae for individual supersymmetric contributions to the 
EDMs due to the Feynman diagrams in Fig.|I[ In our presentation we follow the work 
of Ibrahim and Nath iTll . 









Figure 1: Leading SUSY contributions to the EDMs. The photon and gluon lines are 
to be attached to the loop in all possible ways. 

Neglecting the fiavour mixing, the electromagnetic contributions to the fermion 
EDMs are given by 0: 
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Here 



-.ife 



Mi 



(18) 

with 03 being the gluino phase and Dq defined by D^^M-Dq — dia,g{M~^, M~J. The 
sfermion mass matrix Mj is given by 

Ml^ + m/2 + (I - Qf sin^ 9w) cos 2/3 mf{A} - ^iRf) 

mf{Af - fi*Rf) M^ + nif'^ + Af^Qfsm^ew cos 2/? 

where Rf — cot/9 (tan/?) for = 1/2 (—1/2). The chargino vertex F^jfe is defined as 

r^ife = KaU*^D^,k{VlDl^^ - k^V^D:,,) (19) 

and analogously for the electron; here U and V are the unitary matrices diagonalizing 
the chargino mass matrix: U*M^+V~^ 
Yukawa couplings 



diag(m^+, m^+). The quantities Kf are the 



md,e 



V 2mw sm/3 ' v 2mw^ cos /3 
The neutralino vertex rjfik is given by 

-v^{tan^H^(Q/ - 73,)^!^ + l3,X2i}D},^ - KfXuD%, 
V2taii6wQfXiiDf2k — HfX^iDfH; 



(20) 



Vfik 



X 



(21) 

where I3 is the third component of the isospin, 6 = 3 (4) for 73 = —1/2 (1/2), 
and X is the unitary matrix diagonalizing the neutralino mass matrix: X'^M^oX = 
diag(m^o, m^o, m^o, m^o). In our convention the mass matrix eigenvalues arc positive 
and ordered as m^o > m^o > ... (this holds for all mass matrices in the paper). The 
loop functions A{r), B{r), and C(r) are defined by 



A{r) = 
B{r) = 
C{r) 
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The chromoelectric contributions to the quark EDMs are given by 
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(23) 



Finally, the contribution to the Weinberg operator from the two-loop gluino- 
top-stop and gluino-bottom-sbottom diagrams reads 



d' 



G 



-?>asmt(^\ Im(rf ) ^ B.{zi,Z2,Zt) + {t b) 
KAn / mi 



(24) 



where 



Mr 



where Zi = I ) , zt = i^j . The two-loop function H{zi, Z2, zt) is given by [T^ 



1 N1N2 

H{zi,Z2,Zt) = - dx du dyx{l-x)u 
I Jo Jo Jo 



Ni = u{l — x) + Ztx{l — x){l — u) — 2ux[ziy + 2:2(1 — y)] 



(25) 



'l-x?(l-u? + 
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-x^l-u) 



D = u{l — x) + Ztx{l — x){\ — u) + ux[ziy + ^2(1 — y)] ■ 



(26) 



The numerical behaviour of this function was studied in |T5[. We emphasize that the 
b-quark contribution is significant and often exceeds the top one. 

Before we proceed to the discussion of the EDM suppression mechanisms, let us 
consider the effect of other potentially nonnegligible two-loop contributions. 



3.2 Barr-Zee type EDM contributions 

In view of considerable recent interest in the subject we will consider the two-loop 
Barr-Zee type contributions separately. We will follow the work of Chang, Keung, and 
Pilaftsis 0. 



a 



a 




Figure 2: Barr-Zee type contributions to the EDMs. 



In Ref . ||1 7|| Barr and Zee have presented two-loop Higgs- mediated EDM contribu- 
tions which can be competetive with the Weinberg three-gluon operator. A supersym- 
metric version of the Barr-Zee graphs (Fig.|3) was studied in |16|. In what follows we 
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will analyze only the leading contributions to the EDMs presented in |T^. The EDMs 
arise due to CP- violating couplings of the (s)fermions to the CP-odd Higgs boson oq. 



The EDM and the CEDM of a light fermion / are computed to be [IIB] 



df/e = Qf 
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where i?/ = cot/? (tan/?) for Is = 1/2 (—1/2) and the two-loop function F( 
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The CP-violating couplings ^t,h are given by 
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(29) 



with dii being the standard stop and sbottom mixing angles; bq = Arg(y4g — Rqfi*) and 
v=174 GeV. Note the difference from |]T6| in the definitions of fi and v, also see the 
corresponding Erratum. 

In our numerical analysis, besides assuming the radiative electroweak symmetry 
breaking, we use the following tree level mass of the CP-odd Higgs |]ir 



M 



m 



Hi 



+ m 



H2 



2i/ir, 



(30) 



which is a function of tan f] and other GUT scale input parameters. Strictly speaking, 
this formula is valid for a CP-conserving case, however the EDMs are not very sensitive 
to the exact value of Ml and an inclusion of loop corrections and CP-phases does not 
alter our results. 

In Fig.^ we present a typical Barr-Zee type EDM behaviour as a function of tan/5 
for 0^ = 7r/2, 0^ = 0. The other parameters are fixed to be mo = mi/2 = A = 200 
GeV. The values of tan (3 beyond 42 are not displayed for this parameter set since the 
CP-odd scalar becomes massless in this region and the pattern of the EW symmetry 
breaking becomes unacceptable. We observe that generally these EDM contributions 
by themselves do not impose significant constraints on the GUT scale A-term phases 
even at large tan/5; as can be seen from the plot, the Barr-Zee contributions typically 
are one-two orders of magnitude below the experimental limit. One of the reasons for 
this is that the third generation A-term phases reduce by an order of magnitude due to 
the RG evolution at large tan/?. Also the value of the yU parameter is typically below 
500 GeV owing to the radiative electroweak symmetry breaking. Other factors which 
distinguish our results from those of |T^ are the imposition of Eq. (|30| ) and utilization of 
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the chiral quark neutron model[| For comparison, in Fig.^ we provide the contribution 
of the Weinberg operator which also arises at the two-loop level. 

The constraints become more restrictive for larger A-terms (~ 3mo) and larger 
mo/m3 ratios. In Fig.H we set mo = 500 GeV, mi/2 = 200 GeV, and A = 600 GeV. In 
this case the A-term phase is not "diluted" as much as before and for some parameters 
the Barr-Zee EDMs can be close to the experimental limit. A similar effect can be 
achieved in models with non-universal gaugino masses by introducing a 0{1) gluino 
phase. With the parameter set of Fig.|] the CP-odd Higgs becomes unacceptably light 
around tan P ~ 36. 

/^From the point of view of low energy theory, the Barr-Zee type contributions 
can provide useful constraints on the phases of the third generation A-terms ||16||,|1T9|. 
One can imagine a situation in which the first two generation CP-violating effects are 
suppressed (as in the decoupling scenario), then the EDMs would constrain the third 
generation phases. We find however that typically the Weinberg three-gluon operator 
is considerably more sensitive to such phases and provides more severe constraints even 
at large tan For example, for the parameters of Figs.^ and ^ the contribution of 
the Weinberg operator exceeds that of the Barr-Zee type graphs by one-two orders of 
magnitude. 



Finally, the Z- and W-mediated Barr-Zee type graphs have been analyzed in ^0 
and found to be significantly smaller than those considered above. A number of other 
subleading two loop contributions such as the gluino CEDM-induced quark EDMs, etc. 
have been studied in pi 



Although taken into account, this entire class of diagrams is numerically unimpor- 
tant in our analysis. 



4 Suppression of the EDMs in SUSY models 
4.1 Small CP-phases 

For a light (below 1 TeV) supersymmetric spectrum, the SUSY CP phases have to 
be small in order to satisfy the experimental EDM bounds (unless EDM cancellations 
occur). In Figs.^^ we illustrate the EDMs behaviour as a function of the CP-phases 
in the mSUGRA-type models, where we have set mo = mi/2 = A = 200 GeV. At low 
tan /3, the EDM constraints impose the following bounds (at the GUT scale) : 

< 10"^ - 10"^ , 
0,. < 10-=^ - 10^2 ^ 

<Pgau9. < 10-2 . (31) 

For tan/? > 3, these bounds become even stricter (for 0^ and 4>gaug. the bounds are 
roughly inversely proportional to tan/3). We note that 0a is less constrained than 0^ 

^ We observe similar behaviour in the parton quark model. 

^The Weinberg operator contribution can be suppressed by increasing the gluino mass. However, in 
mSUGRA the Barr-Zee type contributions will also get a suppression factor due to the RG "dilution" 
of the phases of the third generation A-terms. 



10 



and (j)gaug.- There are two reasons for that: first, (j)A is reduced by the RG running from 
the GUT scale down to the electroweak scale and, second, the phase of the {Sfi)LR 
mass insertion which gives the dominant contribution to the EDMs is more sensitive 
to 0^ and 4>gaug. due to \A\ < yLitan/5. 

SUSY models with small CP-phases can be motivated by the approximate CP 
symmetry The well established experimental results exhibit small degree of CP 
violation. Thus it is conceivable that all existing CP-phases are small, including the 
CKM phase. The smallness of CP violation can be explained, for example, by a small 
ratio of the scale at which the CP symmetry gets broken spontaneously and the scale 
at which it is communicated to the observable sector. 

Currently the CKM phase is consistent with zero and it could be supersymmetry 
that is responsible for the observed values of e and e'. This does not require large 
supersymmetric phases. In fact in models with non-universal A-terms e and e' can 
even be saturated with C(10~^) phase of the mass insertion {8^2) lr O- In this case, 
\{5i2) lr\ is required to be (9(10~^) which naturally appears in models with matrix- 
factorizable A-terms of the form B -Ya - C ^ where B and C are flavour matrices. The 
EDM bounds serve to constrain the flavour structure of B and C. Another possibility 
to produce the desired (5i2)l_r is to use asymmetric A-term textures in string-motivated 
models (where the standard supergravity relation Aij = AijYij is assumed). 

Encouragingly, the phase of the /x-term of order 10~^ may be sufficient to produce 



the observed baryon asymmetry ||2J], which is in marginal agreement with the EDM 



bounds. The hypothesis of the approximate CP symmetry is currently being tested in 
the B physics experiments where the Standard Model predicts large CP-asymmetries. 
It is noteworthy that the smallness of the CP-phases in this picture does not constitute 
fine-tuning according to the t'Hooft's criterion p5[ since setting them to zero would 
increase the symmetry of the theory. 

We remark that small CP-phases may also arise due to the dynamics of the system. 
For instance, in weakly coupled heterotic string models, small soft and CKM phases 
arise when the T-moduli get VEV's close to the edge of the fundamental domain 
which is often the case E^. This mechanism however relies on the assumption that 



the dilaton has a real VEV, so this model as it stands does not solve the SUSY CP 
problem. Nevertheless, it may serve as a step toward a consistent string model with 
naturally small CP phases. 

Note that EDMs constrain only the physical phases ([T5|). One can imagine a 
situation when the individual phases are 0{1) whereas the physical ones are small. This 
occurs for example in gauge mediated SUSY breaking models, see |^ and references 
therein. 



4.2 Heavy SUSY scalars 

This possibility is based on the decoupling of heavy supersymmetric particles. Even if 
one allows 0{1) CP violating phases, their effect will be negligible if the SUSY spectrum 
is sufficiently heavy . Generally, SUSY fermions are required to be lighter than the 
SUSY scalars by, for example, cosmological considerations. So the decoupling scenario 
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can be implemented with heavy sfermions only. Here the SUSY contributions to the 
EDMs are suppressed even with maximal SUSY phases because the squarks in the loop 
are very heavy and the mixing angles are small. 

In Fig.^we display the EDMs as functions of the universal scalar mass parameter mo 
for the mSUGRA model with maximal CP-phases 0^ = 0a = 7r/2 and mi/2 = A = 200 
GeV. We observe that all EDM constraints except for that of the electron require tuq 
to be around 5 TeV or more. The mercury constraint is the strongest one and requires 

decoupl . 

~ 10 TeV . (32) 

This leads to a serious fine-tuning problem. Recall that one of the primary motivations 
for supersymmetry was a solution to the naturalness problem. Certainly this motivation 
will be entirely lost if a SUSY model reintroduces the same problem in a different sector, 
i.e. for example a large hierarchy between the scalar mass and the electroweak scale. 

The degree of fine-tuning can be quantified as follows. The Z boson mass is deter- 
mined at tree level by 

1 _ mjj^-ml^tan-'P . 

2^^- tan2 /3 - 1 ^ ■ ^^^^ 



One can define the sensitivity coefficients [28 



C; = 



9 In m| 



d In a? 



(34) 



where are the high energy scale input parameters such as vrixji^ ""^O) etc. Note that 
[I is treated here as an independent input parameter. A value of Cj much greater than 
one would indicate a large degree of fine-tuning. The Higgs mass parameters are quite 
sensitive to mo, so for mo=10 TeV we find 

Cmo ~ 5000 (35) 

for the parameters of Fig.^. For the universal scalar mass of 5 (3) TeV the sensitivity 
coefficient reduces to 1300 (500). This clearly indicates an unacceptable degree of 
fine-tuning. 

This problem can be mitigated in models with focus point supersymmetry, i.e. when 
m^jj is insensitive to mo I^Qj. However, this mechanism works for mo no greater than 



2-3 TeV which is not sufficient to suppress the EDMs. Another interesting possibility is 
presented by models with a radiatively driven inverted mass hierarchy, i.e. the models 
in which a large hierarchy between the Higgs and the first two generations scalar masses 
is created radiatively However, a successful implementation of this idea is far from 
trivial One can also break the scalar mass universality at the high energy scale 
321. In this case, either a mass hierarchy appears already in the soft breaking terms 



or certain relations among the soft parameters must be imposed (for a review see |]33|). 
These significant complications disfavour the decoupling scenario as a way to solve the 
SUSY CP problem, yet it remains a possibility. 

Note that the decoupling scenario rules out supersymmetry as a possible explana- 
tion of the recently observed 2.6 cr deviation in the muon g — 2 from the SM prediction 
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34| . This scenario may also lead to cosmological difficulties, in particular with the 
relic abundance of the LSPs since the LSP annihilation cross section falls rapidly as 
nisfermion iucreases. Concerning the other phenomenological consequences, we remark 
that the SUSY contributions to the CP-observables involving the ffist two generations 
(such as e, e') are negligible, although those involving the third generation may be 
considerable. The corresponding CP-phases are constrained through the Weinberg op- 
erator contribution to the neutron EDM, which typically prohibits the maximal phase 



= 7r/2 while still allowing for smaller Oil) phases (Fig| 

GUT 



4.3 EDM cancellations 

The cancellation scenario is based on the fact that large cancellations among differ- 
ent contributions to the EDMs are possible in certain regions of the parameter space 
[0 , , , which allow for 0{1) flavour — independent CP phases. Although not 



particularly well motivated, this possibility is interesting since, when supplemented 
with a real non-trivial flavour structure, it allows for significant supersymmetric con- 
tributions to CP observables including those in the B system. In this case, the SM 
predictions can be significantly altered leading for instance to a nonclosure of the uni- 
tarity triangle (see M. Brhlik et al. in ||23|). Given the appropriate A-terms' flavour 
structures, the parameters e and e' can be of completely supersymmetric nature 



Large flavour-independent SUSY phases may also be responsible for electroweak baryo- 
genesis Thus the cancellation scenario presents a interesting alternative to the 

decoupling and approximate CP solutions. 

For the case of the electron, the EDM cancellations occur between the chargino 
and the neutralino contributions. For the case of the neutron and mercury, there are 
cancellations between the gluino and the chargino contributions as well as cancellations 
among contributions of different quarks to the total EDM. A number of approximate 
formulae quantifying the cancellations are presented in Refs.[^,p6[]. 



In what follows we examine the cancellation scenario in the universal and nonuni- 
versal cases, with the latter being motivated by Type I string models. We note that 
the CP phases are to be understood modulo it. 



4.3.1 EDM cancellations in mSUGRA-type models. 

These are mSUGRA-type models allowing different phases for different gaugino masses 
while keeping their magnitudes the same. As we will see, mSUGRA models with zero 
gaugino phases can realize the cancellation scenario. An introduction of the gaugino 
phases makes the cancellations much harder to achieve (if possible at all) and leads to 
a much higher degree of fine-tuning. The parameters allowing the EDM cancellations 
strongly depend on the neutron model. For example, in the parton model, it is more 
difficult to achieve these cancellations due to the large strange quark contribution. 
Therefore, one cannot restrict the parameter space in a model-independent way and 
caution is needed when dealing with the parameters allowed by the cancellations. 
In mSUGRA, the EDM cancellations can occur simultaneously for the electron. 
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neutron, and mercury along a band in the (0a, 0/^) plane ( Figs. |T0| and |TT]) . However, in 
this case the mercury constraint requires the /x phase to be (9(10~^) and the magnitude 
of the A-terms to be suppressed (~ O.lmo) which results in only a small effect of the 
A-terms on the phase of the corresponding mass insertion. This is to be contrasted 
with simultaneous EEDM/NEDM cancellations which allow for 0^ ~ (9(10^^) and 
A ~ mo In that case the individual contributions to the EDMs exceeded the 



experimental limit by an order of magnitude (or more). Owing to the addition of the 
mercury constraint, the EDM cancellations become much milder as illustrated in Fig.|l2|. 
For example, without these cancellations the EDMs would exceed the experimental 
limit only by a factor of a few. Obviously, the border between the cancellation and the 
small phases scenarios becomes blurred. This is even more so at larger tan/3 (Fig.^), 
in which case the cancellation band becomes narrower and the limits on 0^ become 
tighter. We note that there could exist some points which allow for (p^ ~ O{10~^] 
however these are rare exceptions and such points do not form a band. 



As noted in Ref.|]lT|, in the case of zero gaugino phases, the band in the 



plane where the cancellations occur can approximately be described by a relation 



V - 



— asin0A , (36) 



where a is a constant depending on the parameters of the model which represents 
the maximal allowed phase 0^. For example, for the chiral quark neutron model and 
tan/? = 3, mo = mi/2 = 200 GeV, and \A\ = 40 GeV (parameters of Fig.[l|), a = 0.017. 
Of course, the value of a depends on the input parameters such as the quark masses, 
the GUT scale, the SUSY scale, etc. and also involves numerical uncertainties, so 
caution is needed when treating this number. The maximal allowed 0^ is roughly 
inversely proportional to tan/5, e.g. for tan/3 = 10 we have a = 0.005 (Fig.|l^). In 
the case of the parton model, the cancellations occur, for example, with a = 0.006 and 
tan/3 = 3, mo = mi/2 = 200 GeV, and \A\ = 20 GeV (Fig.|ll|). As mentioned earlier, 
the cancellations are harder to achieve in the parton model, which results in tighter 
limits on 0^ (~ O{10^^)). In both cases 0^ is restricted to be of order 10^'^ — 10~^, 
whereas 0a can be arbitrary (however physical effects due to 0a will be suppressed 
because of the small magnitude of A). The GUT parameters given above imply that 
the squark masses at low energies are about 500 GeV, so these bounds can be relaxed 
only if the squark masses are over 1 TeV. 

If the gluino phase is turned on, simultaneous EEDM, NEDM, and mercury EDM 
cancellations are not possible. The gluino phase affects the NEDM cancellation band 
by altering the relation ([36|): 



0^ ~ -a sin(0A + a) - c , (37) 

while leaving the EEDM cancellation band almost intact (Fig.|14D. An introduction 
of the bino phase 0i qualitatively has the same "off-setting" effect on the EEDM 
cancellation band as the gluino phase does on that of the NEDM (Fig.|T5|). Note that 
the bino phase has no significant effect on the neutron and mercury cancellation bands 
since the neutralino contribution in both cases is small. When both the gluino and 
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bino phases are present (and fixed), simultaneous electron, neutron, and mercury EDMs 
cancellations do not appear to be possible along a band. The reason is that the mercury 
EDM depends on 03 more strongly than the NEDM does, so that an introduction of 
the gluino phase splits the bands where the mercury and neutron EDM cancellations 
occur as illustrated in Figs.O and IIR Thus we see that zero gaugino phases are much 



more preferred by the cancellation scenario. 

The cancellation scenario involves a significant fine-tuning. Indeed, restricting the 
phases to the band where the cancellations occur does not increase the symmetry of the 
model and thus is unnatural according to the t'Hooft's criterion ||2^. It is a non-trivial 
task to quantify the degree of fine-tuning in this case. One possibility is to define an 
EDM sensitivity coefficient with respect to the CP-phases, in analogy with Eq.^. We 
typically find it to be between 30 and 100 on the cancellation band. This however 
represents only a local behaviour of the EDMs. In other words it shows how easy it is 
to spoil the cancellations without a reference to how improbable it is to achieve such 
cancellations in the first place. An alternative way to quantify the degree of fine-tuning 
is simply to estimate the probability that a random small area in the (0^, 0^) plane will 
satisfy the cancellation condition. Since any point is not prefered over any other point 
by the underlying theory, this should give a fairly good idea of the minimal degree of 



fine-tuning needed. From Fig.|TO| we obtain 

fine — tuning ~ (probability)"^ > 10^ . (38) 

Note that this estimate does not take into account the fine-tuning of the other soft 
breaking parameters which is necessary to allow for simultaneous EEDM, NEDM, and 
mercury EDM cancellations. Other estimates give a similar number for the universal 
case, whereas for the nonuniversal case the degree of fine-tuning drastically increases: 
only one out of 10^ random points in the parameter space satisfies the cancellation 
condition 



4.3.2 EDM cancellations in D-brane models. 



Let us first briefly review basic ideas of D-brane models (see also Refs. ||39[| and |^). 
Recent studies of type I strings have shown that it is possible to construct a number 
of models with non-universal soft SUSY breaking terms which are phenomenologically 
interesting. Type I models can contain 9-branes, 5j-branes, 7j-branes, and 3-branes 
where the index i = 1,2,3 denotes the complex compact coordinate which is included 
in the 5-brane world volume or which is orthogonal to the 7-brane world volume. 
However, to preserve = 1 supersymmetry in D = A not all of these branes can be 
present simultaneously and we can have (at most) either D9-branes with D5j-branes 
or D3-branes with DTj-branes. 

Gauge symmetry groups are associated with stacks of branes located "on top of 
each other". A stack of A^ branes corresponds to the group U{N). The matter fields 
are associated with open strings which start and end on the branes. These strings 
may be attached to either the same stack of branes or two different sets of branes 
which have overlapping world volumes. The ends of the string carry quantum numbers 
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associated with the symmetry groups of the branes. For example, the quark fields have 
to be attached to the U{3) set of branes, while the quark doublet fields also have to be 
attached to the U{2) set of branes. Given a brane configuration, the Standard Model 
fields are constructed according to their quantum numbers. 

The SM gauge group can be obtained in the context of D-brane scenarios from 
f/(3) X f/(2) X U{1), where the U{3) arises from three coincident branes, U{2) arises 
from two coincident D-branes and U{1) from one D-brane. As explained in detail in 



Ref.|4^, there are different possibilities for embedding the SM gauge groups within 
these D-branes. It was shown that if the SM gauge groups come from the same set 
of D-branes, one cannot produce the correct values for the gauge couplings aj{Mz) 
and the presence of additional matter (doublets and triplets) is necessary to obtain 
the experimental values of the couplings |^ . On the other hand, the assumption that 
the SM gauge groups originate from different sets of D-branes leads in a natural way 
to intermediate values for the string scale Ms — 10^°^^^ GeV EDI. In this case, the 



analysis of the soft terms has been done under the assumption that only the dilaton 
and moduli fields contribute to supersymmetry breaking and it has been found that 
these soft terms are generically non-universal. The MSSM fields arising from open 
strings are shown in Fig.|^. For example, the up quark singlets u'^ are states of the type 
C^^^, the quark doublets are C^^^ , etc. The presence of extra {Dq) branes which are 
not associated with the SM gauge groups is often necessary to reproduce the correct 
hypercharge and to cancel non- vanishing tadpoles. 




Figure 3: Embedding the SM gauge group within different sets of D-branes. The extra 
Dg brane (52) is marked by a cross. 

Recently there has been a considerable interest in supersymmetric models derived 



from D-branes |^,[^]. In a toy model of Ref . P2| , the gauge group SU{3)c x f/(l)y 
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was associated with 5i branes and SU{2)l was associated with 52 branes. It was 
shown that in this model the gaugino masses are non-universal (Mi = M3 7^ M2) so 
that the physical CP phases are (f)i = (f>3, 4>a and 0^. Non- universality of the gaugino 
masses allowed to enlarge the regions of the parameter space where the NEDM/EEDM 
cancellations occured. However, such a model is oversimplified and the relation 0i = 
03 7^ 02 does not appear to hold in more realistic models In what follows we will 



consider the EDM cancellations in this and a more realistic D-brane models. 

The model in which f/(3), U{2), and U{1) originate from different sets of branes 
is phenomenologically interesting. In this case one naturally obtains an intermediate 
string scale (10^° — 10^^ GeV), although higher values up to 10^^ GeV are still al- 
lowed. Both the up and the down type Yukawa interactions are allowed, while that 
for the leptons typically vanishes (depending on further details of the model) P0[ . 
The (anomaly-free) hypercharge is expressed in terms of the U{1) charges Qi,2,3 of the 
C/(l)i,2,3 groups: 

Y = -Iq3-^Q2 + Qi, (39) 
with the following (^3,^27^1) charge assignment: 

q = (1, -1, 0) , = (-1, 0, -1) , = (-1, 0, 0) , (40) 
/ = (0,1,0), e^ = (0,0,1), 
i/2 = (0,1,1), ifi = (0,1,0). 

The gaugino masses in this model are given by 

M3 = V3m3/2sin^ e-*"% (41) 



M2 = v^m3/2 ©1 cos6' e" 
My = v^m3/2 ayiMs) 



2 12 

-63 COS ^e"'°3 + — -61 cos 0e~'"^ + sin ( 



ai(Ms) a2{Ms) 3^3 (M5) 



where 



1 2 
+ ^TTT + T^^TTT • (42) 



ay^Ms) ai{Ms) a2{Ms) ?>a3{Ms) 



Here correspond to the gauge couplings of the U{k) branes. As shown in Ref. [40 



ai{Ms) ~ 0.1 leads to the string scale Ms ~ 10^^ GeV. The parameters 9 and Gj 



parameterize supersymmetry breaking in the usual way [p9| : 

= V3(S + 5*)m3/2 sin^ e"*"^ , 
r = V3{Ti + T*)ms/2 cos 9 e^e"'"^ , (43) 

and i = 1,2,3 labels the three complex compact dimensions. 
The soft scalar masses are given by 



2 2 

= ^3/2 
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(44) 



and the trilinear parameters are 



A-u 


V3 

= ^"^3/2 


"(026-^"^ 


- Qie-'"' 


- 036"*"='^ 


COS 6 




A, 


= ^"^3/2 


"(036-'"^ 


- Q^e-'"' 


- 026"^"^^ 


COS 9 


- sin 9 e-^"= 


Ae 


= . 













(45) 

(46) 
(47) 

We observe that the angles 0j and 9 are quite constrained if we are to avoid negative 
mass-squared 's for squarks and sleptons. In what follows we set 03 = and ai = 02; 
then the soft terms are parameterized in terms of the phase = ai — a^. 

In Fig.|l^ we display the bands allowed by the electron (red), neutron (green), 



and mercury (blue) EDMs. In this figure, we set m 



3/2 



150 GeV, tan/? 



92 = 92 = 1/2^ cos^^ = 2sin2^ = 2/3, and ai{Ms) ~ 1 with Ms being the GUT 
scale. For the plot to be more illustrative, we do not impose any additional constraints 
besides the EDM ones (i.e. bounds on the chargino and slepton masses, etc.). It is 
clear that even though simultaneous EEDM/NEDM cancellations allow the phase 
to be 0{1), an addition of the mercury constraint requires all phases to be very small 
(modulo tt) and thus practically rules out the cancellation scenario in this context. 
The situation becomes even worse in the case of an intermediate string scale ~ 10^^ 
GeV (i.e. ai{Ms) ~ 0.1), see Fig.|l^. We find it quite generic that the mercury EDM 
behaviour in D-brane models is very different from that of the electron and neutron 
and thus is crucial in constraining the parameter space. The major difference from 
the mSUGRA-type models with fixed 0y and 03 is that the phase of the A-terms is 
correlated with the gaugino phases resulting in the cancellation bands which are not 



described by a simple relation 



—a sim (f)A + a) — c. 



Next we consider the model of Ref . ||4^ . The (corrected) soft terms for this model 
read (for 93 = 0) 

My = Ms = -A = V3ms/2 cos9 016"^°^ , 
M2 = y3m3/2Cos^02e-'°2 , 



m, 



m 



3/2 



m 



R 



"^3/2 ( 



;i-^sin^^); 
1 -Scos^e 0^) 



(4^ 
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To illustrate the EDM constraints, we choose the parameters which allow for simulta- 
neous EEDM/NEDM cancellations, namely m^i2 = 150 GeV, tan/5 = 2, Oi = 0.9, and 



6 = 0.4 as given in Ref.|^2[. Fig.|T8| shows that the mercury constraint has the same 
behaviour as in the model considered above and rules out large CP-phases. 

We see that the cancellation scenario in simple models faces a number of difficulties. 
Presently available string-motivated models with non-universal gaugino masses cannot 
accommodate simultaneous electron, neutron, and mercury EDM cancellations. In the 
mSUGRA, such cancellations are possible but require a significant fine-tuning. The 
addition of the mercury EDM bound restricts the phase of the yU term to be O{10~^) if 
we are to achieve the EDM cancellations along a band in the {(pA, <P^i) plane. Without 
an additional SUSY flavour structure, the CP-phases allowed by the cancellations will 
have very small observable effects. Even in a more general situation (unconstrained 
MSSM), the phases allowed by the EDM cancellations typically lead to small CP- 
asymmetries ( < 1%) in collider experiments |3^. Testing the cancellation scenario 
experimentally may prove to be a challenge. 



4.4 Flavour-off-diagonal CP violation 

This is one of the more attractive possibilities to avoid overproduction of EDMs in 
SUSY models. Nonobservation of EDMs may imply that CP-violation has a flavour- 
off-diagonal character just like in the Standard Model. The origin of CP- violation 
in this case is closely related to the origin of the flavour structures rather than the 
origin of supersymmetry breaking. While models with flavour-off- diagonal CP violation 
naturally avoid the EDM problem, they have testable effects in K and B physics. 

This class of models requires hermitian Yukawa matrices and A-terms, which forces 
the flavour-diagonal phases to vanish (up to small RG corrections) in any basis. The 
flavour-independent quantities such as the /x-term, gaugino masses, etc. are real. This is 
naturally implemented in left-right symmetric models and models with a horizontal 
flavour symmetry |^ . 



In the left-right models, the hermiticity of the Yukawas and A-terms as well as the 
reality of the /x-term is forced by the left-right symmetry. The SU{2)l gaugino mass 
is in general complex, so in order to suppress the EDMs the additional assumption 
of its reality is needed. The phenomenology of such models in the context of up- 
down uniflcation has been studied in |^6[. The left-right symmetry appears to be too 



restrictive in this case to satisfy all of the phenomenological constraints; however the 
decisive test will be undertaken at the B factories. 

Another possibility is based on a horizontal U{3)h sjTiimetry Hermitian 
Yukawa matrices may appear due to a (gauged) horizontal symmetry U{3)h which gets 
broken spontaneously by the VEVs of the real adjoint flelds (a = 1..9; a = u,d, I, ..) 
Some of these real VEVs also break CP since some of the components of are 



CP-odd. As a result, CP violation appears in the superpotential through complex 
Yukawa couplings, whereas the /i-term is real since it arises from a U{3)h invariant 
combination of the type T^T^. An effective f/(3)i^-invariant superpotential of the type 
W = ^HiQi{T^\%jbj produces the Yukawa matrix {Y%j = 3^{T^){\%j, where 
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A° is proportional to the unit matrix and A^~^ are the Gell-Mann matrices. Note that 
the real fields may only come from a non-supersymmetric (anti-brane) sector. If 

are the scalar components of chiral multiplets, they are intrinsically complex and 
hermitian Yukawas in this case arise only if the VEVs of are real. 

The gaugino masses in this model are not forced to be real by symmetry. One has 
to make an additional assumption that either the gaugino masses are universal and 
the corresponding phase can be rotated away or that the SUSY breaking dynamics 
conserve CP which seems natural if CP breaking is associated with the origin of fiavour 
structures. In other words, the SUSY breaking auxiliary fields get real VEVs as a result 
of the underlying dynamics such as the dilaton stabilization in Type I string models 
|9| or the effective potential minimization in heterotic string models [5C]. Further, 
if the Kahler potential is either generation independent or left-right symmetric, the 
A-terms are hermitian This leads to very small phases in the fiavour-diagonal 

mass insertions which are responsible for generating EDMs and the EDM bounds are 
easily satisfied. 

Both the left-right model and the model with a horizontal symmetry mitigate the 
strong CP-problem (under the additional assumptions given above). The 6 parameter 
vanishes at both the tree and the leading log one-loop levels. However neither of these 
models solves the problem completely due to the significant one loop finite corrections 
which appear mostly due to a non-degeneracy of the squark masses and a misalignment 



between the Yukawas and the A-terms 51 



In this class of models we have the following setting at the high energy scale: 



Arg(Mfe) = Arg(/i) = . 



(49) 



Generally, the off-diagonal elements of the A-terms can have 0{1) phases without 
violating the EDM constraints. Due to the RG effects, large phases in the soft trilinear 
couplings involving the third generation generate small phases in the flavour-diagonal 
mass insertions for the light generations, and thus induce the EDMs. For example, the 
A-terms of the formQ 



( 1 

a 



12 



ai2 
1 
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a 



23 



ai3^ 

^23 
^33/ 



and the following GUT-scale hermitian Yukawa matrices 



(50) 



( 4.1x10-^ 6.9 X 10-M -1.4x10-2 

-1.4 X 10"^ i 



-6.9 X 10-* i 3.5 X 10" 



V 



-1.4x10-2 1.4xlO-M 6.9x10-1 



/ 1.3 X 10-^ 

2.0 X 10-^ - 1.8 X 10-^ 
V -4.4 X 10-^ 



2.0 X 10-^ + 1.8 X 10-^ i -4.4 x 10"^ ~ 



9.3 X 10-^ 
-7.0 X 10-^ i 



7.0 X 10-* i 
1.9 X 10 



-2 



•^The standard supergravity relation (^Aa)ij 
sumed. 



i,J^a)ij{ya)ij for the trilinear soft couplings is as- 
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typically induce the mercury EDM of the order of the experimental limit if 033 is of 
order 1, whereas the induced NEDM is 1-2 orders of magnitude below the experimental 
limit. If one uses Yukawa textures with smaller Y-^, this RG effect will be suppressed 
rendering the induced mercury EDM far below the experimental limit. 

The supersymmetric contribution to the e' parameter is suppressed in models with 
hermitian flavour structure^ [^. This occurs due to severe cancellations between 
the contributions involving {Sf2) lr and (5i2)rl niass insertions (we use the standard 
definitions of |^). Due to the hermiticity (5i2)Li? — {^12) RL, whereas they contribute 



Im 



< 10~^ which 
On the other 



to e' with opposite signs. Typically we find 
produces e' an order of magnitude below the experimental limit 
hand, similar cancellations do not occur for the e parameter and the SUSY contribution 
to e can be even dominant. The value of (5^2)^^? which saturates the observed |£:| ~ 



2.26 X 10-'^ is given by y |Im((5f2)iij| — "^-^ x 10 ^ for the gluino and squark masses of 
500 GeV g. 

For completeness, we provide the bounds on the imaginary parts of the mass inser- 
tions 

i5f^)LR = ^Mic^Mhvm - >^/^^V^2(l)) (51) 

fit 

derived from the leading gluino contributions to the electromagnetic (NEDM) operator 
and the bino contribution to the electron EDM. We update the bounds of Ref. and 
also include the QCD correction factor. The advantage of the mass insertion approach 
is that it allows to obtain model independent bounds and thus is quite useful when 
dealing with complicated flavour structures. For comparison we present the bounds for 
the chiral (Table 1) and the parton (Table 2) neutron models. 

A drastic improvement of the bounds comes from the addition of the mercury 
EDM constraint. Using the expressions for the chromomagnetic moments of Ref. |p7||. 



in Table 3 we present the bounds on the mass insertions from the gluino contributions to 
the mercury EDM. For \lm{Sfi)Lf{)\ and |Im((5"]^)i^)| these bounds turn out to be very 
strict, more than an order of magnitude stricter than those imposed by the NEDM. This 
severely restricts the CP-asymmetry Acp{b — >■ 57) in models with hermitian flavour 
structures. The reason is that in order to obtain a large (~ 10%) SUSY contribution 
to this observable, the elements of the A-terms involving the third generation have to 
be larger than 1. This induces via the RG running a considerable \lm{6f'^)LR)\, often 
in conflict with the bounds of Table 3. We find that with the above Yukawa textures 
Acp{b —>■ 57) is allowed to be no more than 2-3%. One can relax this constraint by 
using different hermitian textures, especially with suppressed Y{^. In this case the 
CP-asymmetry can be as large as 6-7%. 

To conclude this section, we remark that the problem of baryogenesis in this class 
of models requires careful investigation and at the moment it is unclear whether or not 
large flavour off-diagonal SUSY phases can produce sufficient baryon asymmetry of the 
universe. 



"^For phenomenology of models with non-universal (but non-hermitian) A-terms see also 
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5 Discussion and conclusions 



We have systematically analyzed constraints on supersymmetric models imposed by the 
experimental bounds on the electron, neutron, and mercury electric dipole moments. 
We find that the EDMs can be suppressed in SUSY models with 

1) small SUSY CP-phases ( < 10"^). This possibility can be motivated by the 
approximate CP-symmetry which also implies that the CKM phase is small. This 
provides testable signatures for B-factories' experiments. 

2) heavy SUSY scalars, 

fyisfermion ~ 10 TeV. In this class of models there is a large 
hierarchy between the SUSY and electroweak scales, which is hard to realize without 
an extreme fine-tuning. 

3) EDM cancellations. We have analyzed the possibility of such cancellations in 
D-brane and mSUGRA-like models with nontrivial gaugino phases. We find that, with 
the addition of the mercury EDM constraint, only the EDM cancellations in mSUGRA 
(01 ~ 03 ~ 0) survive in any considerable part of the parameter space. Even in this case 
the cancellations require small 0^ ~ 10~^ and suppressed |y4| (~ O.lmo). As a result, 
the border between the small phases and the cancellation scenarios fades away. In 
addition to the finetuning problem, models with the EDM cancellations lack predictive 
power as it is unclear whether the allowed CP-phases can have observable effects. 

4) fiavour-off-diagonal CP violation. This can occur in models with CP-conserving 
SUSY breaking dynamics and hermitian flavour structures. Such models allow for 0{1) 
flavour off-diagonal phases which can have signiflcant effects in K and B physics. 

It is also possible to combine different mechanisms to suppress the EDMs. For 
example, a "hybrid" of the decoupling and the cancellation scenarios was considered 



in Ref.fSJ]. Such models seem to share shortcomings of both "parents" without an 
apparent advantage over either of them. 

There is also a rather radical proposal that all gaugino masses and the A-terms are 
vanishingly small |j55| . Clearly this eliminates all of the physical phases in Eq. ([I5|) and 
thus produces no CP violation. Such a strong assumption requires a flrm motivation 
such as the continuous R-symmetry. However, by the same token, the R-symmetry 
eliminates the S/i term thereby leading to a very light CP-odd Higgs boson. Such a 
scenario faces a number of difficulties with experimental results. 

To summarize, we have studied the supersymmetric CP problem taking into ac- 
count all of the current EDM constraints. Our conclusion is that there remain two 
attractive ways to avoid overproduction of the EDMs in SUSY models. The flrst pos- 
sibility is that CP is an approximate symmetry of nature and the second one is that 
CP violation has a flavour off-diagonal character just as in the Standard Model. 
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Tabic 1: Bounds on the imaginary parts of the mass insertions. The chiral quark 
model for the neutron is assumed. Here x = m'^/m'^ = rnj^/rrij, rriq = 500 GeVjUij = 
100 GeV. For different squark/slepton masses the bounds are to be multiplied by 
mg/500 GeV or m^/lOO GeV. 
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Table 2: Bounds on the imaginary parts of the mass insertions for the parton neutron 
model. For the squark masses different from 500 GeV, the bounds are to be multiphed 
by mg/500 GeV. 
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Table 3: Bounds on the imaginary parts of the mass insertions imposed by the mercury 
EDM. For the squark masses different from 500 GeV, the bounds are to be multiphed 
by mg/500 GeV. 
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Figure 4: Barr-Zee and Weinberg operator induced EDMs as a function of tan/3. 1 
- electron Barr-Zee EDM, 2 - neutron Barr-Zee EDM, 3 - neutron EDM due to the 
Weinberg operator. Here uiq = mi/2 = A = 200 GeV, 0^ = 0, 0^ = 7r/2. 
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Figure 5: Barr-Zee and Weinberg operator induced EDMs as a function of tan/?. 1 
- electron Barr-Zee EDM, 2 - neutron Barr-Zee EDM, 3 - neutron EDM due to the 
Weinberg operator. Here mo = 500 GeV, mi/2 = 200 GeV, A = 600 GeV, = 0, 
cpA = vr/2. 
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Figure 6: EDMs as a function of (pA- 1 ^ electron, 2 - neutron (chiral model), 3 - 
neutron (parton model), 4 - mercury. The experimental limit is given by the horizontal 
line. Here tan/3 = 3, mo = mi/2 = A = 200 GeV. 
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Figure 7: EDMs as a function of 0^. 1 - electron, 2 - neutron (chiral model), 3 - 
mercury, 4 - neutron (parton model). The experimental limit is given by the horizontal 
line. Here tan/3 = 3, mo = mi/2 = A = 200 GeV. 
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Figure 8: EDMs as a function of the gluino phase 03. 1 - electron, 2 - neutron (chiral 
model), 3 - neutron (parton model), 4 - mercury. The experimental limit is given by 
the horizontal line. Here tan/? = 3, mo = mi/2 = A = 200 GeV. 




Figure 9: EDMs as a function of the universal mass parameter mo. 1 - electron, 2 - 
neutron (chiral model), 3 - neutron (parton model), 4 - mercury. The experimental 
limit is given by the horizontal line. Here tan/? = 3, mi/2 = A = 200 GeV, 0^ = 0a = 
7r/2. 
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Figure 10: Phases allowed by simultaneous electron, neutron, and mercury EDM can- 
cellations in mSUGRA. The chiral quark neutron model is assumed. Here tan/3 = 3, 



mo = mi/2 



200 GeV, A = 40 GeV. 
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Figure 11: Phases allowed by simultaneous electron, neutron, and mercury EDM can- 
cellations in mSUGRA. The parton quark neutron model is assumed. Here tan/? = 3, 
mo = mi/2 = 200 GeV, A = 20 GeV. 
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Figure 12: Chargino - neutralino cancellations for the EEDM. 1 - neutralino, 2 - 



chargino, 3 - total EDM. The SUSY parameters are the same as for Fig.|TO|, which 
allow for simultaneous electron, neutron, and mercury EDM cancellations. 
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Figure 13: Phases allowed by simultaneous electron, neutron (chiral model), and 
mercury EDM cancellations in mSUGRA at larger tan/3. Here tan/5 = 10, mo = 
mi/2 = 200 GeV, A = AO GeV. 
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Figure 14: Bands allowed by the electron (1), neutron (2), and mercury (3) EDMs 
cancellations in the mSUGRA-type model with a nonzero gluino phase. Here tan /3 = 3, 
mo = mi/2 = 200 GeV, A = AO GeV, 03 = tt/IO. 
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Figure 15: Bands allowed by the electron (1), neutron (2), and mercury (3) EDMs 
cancellations in the mSUGRA-type model with nonzero gluino and bino phases. Here 
tan/3 = 3, mo = mi/2 = 200 GeV, A = AO GeV, 0i = 03 = tt/IO. 
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Figure 16: Bands allowed by the electron (1), neutron (2), and mercury (3) EDMs 
in the D-brane model. Here tan/5 = 3, m3/2 = 150 GeV, = &l = 1/2, cos^ 6' = 
2 sin^ e = 2/3, and Ms ~ 10^^ GeV. 




Figure 17: Bands allowed by the electron (1), neutron (2), and mercury (3) EDMs 
in the D-brane model with an intermediate scale. Here tan/3 = 3, 1713/2 = 150 GeV, 
92 = 92 = 1/2, cos2 6 = 2 sin^ 9 = 2/3, and Ms ~ lO^^ Q^y 
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Figure 18: Bands allowed by the electron (1), neutron (2), and mercury (3) EDMs in 
the model of Ref.[||. Here tan/? = 2, mg/a = 150 GeV, 0i = 0.9, 6 = 0.4. 
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